Systems and methods for 3-dimensional interferometric microscopy

ABSTRACT

In one embodiment, an apparatus comprises an optical system with multiple detectors and a processor. The optical system is configured to produce images of an optical source in a first dimension and a second dimension substantially orthogonal to the first dimension at each detector at a given time. Each image from the images is based on an interference of an emission from the optical source in a first direction and an emission from the optical source in a second direction different from the first direction. The processor is configured to calculate a position in a third dimension based on the images. The third dimension is substantially orthogonal to the first dimension and the second dimension.

CROSS-REFERENCE TO RELATED APPLICATION

This application claims priority to U.S. Patent Application Ser. No. 60/871,366, filed on Dec. 21, 2006 and entitled “System and Methods for 3-Dimensional Interferometric Microscopy,” and U.S. Patent Application Ser. No. 60/908,307, filed on Mar. 27, 2007 and entitled “System and Methods for 3-Dimensional Interferometric Microscopy;” both of which are incorporated herein by reference.

BACKGROUND

The invention relates generally to interferometric microscopy including, for example, interferometric microscopy of samples comprised of or labeled with spatially resolvable point sources such as samples that internally radiate light from photo-activated fluorescent protein molecules and/or interferometric microscopy where a phase difference is measured to determine position information in the third dimension (e.g., the z coordinate) while simultaneously measuring position information in the two other dimensions (e.g., the x and y coordinates).

Microscopy systems that measure interferometric phase are known. For example, one known microscopy system measures interferometric phase by illuminating an extended surface of an object with an external light source. Light reflected from the sample interferes with a reference beam in a beam splitter. The interfered beam is imaged onto a detector. The optical path difference between the reference beam and the sample beam can modulate the amplitude of the output interfered beam and provide a measure of an object height.

Also known are microscopy systems that interfere two beams from a sample via opposing objectives to measure only one interfered output beam. Such microscopy systems do not explicitly measure phase. These microscopy systems can resolve structure in the axial dimension.

Such known microscopy systems, however, suffer several disadvantages. For example, such systems cannot be applied to fluorescently-labeled samples that are common in biology because the internally-supplied fluorescent radiation has no useable phase relation with respect to any externally-supplied reference beam used for excitation of fluorescence. Without a reference beam, interference is not achievable. In addition, a single or sequentially phased interferometric measurement cannot be used to separate interferometric phase from amplitude if the radiation is transient or strongly time variant, such as is the case in blinking fluorescent molecules or other single photon sources. Thus, no quantitative measure of phase and positional information can be derived from a single interferometric measurement or a time sequence of measurements. Finally, lateral resolution for far-field interferometric systems is limited by the Abbe diffraction length λ/2NA, where λ is the wavelength and NA is the numerical aperture. Thus, a need exists for an improved microscopy system.

SUMMARY

In one embodiment, an apparatus comprises an optical system with multiple detectors and a processor. The optical system is configured to produce images of an optical source in a first dimension and a second dimension substantially orthogonal to the first dimension at each detector at a given time. Each image from the images is based on an interference of an emission from the optical source in a first direction and an emission from the optical source in a second direction different from the first direction. The processor is configured to calculate a position in a third dimension based on the images. The third dimension is substantially orthogonal to the first dimension and the second dimension.

In some embodiments, molecules of interest that are tagged with fluorescent labels can be localized in 3 dimensions to an accuracy better than the diffraction limit. In such embodiments, the phase information associated with the image can be used to determine information about the position of a molecule in the z coordinate such that a 3D (three-dimensional) representation with sub-diffractive accuracy <<λ/2 in the 3 dimensions can be achieved. In other words, information about the position of a molecule in the third dimension (e.g., the z coordinate) can be obtained while information about the position of the molecule in the two other dimensions (e.g., the x and y coordinates) can be obtained.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a system block diagram of a two-phase interferometric microscopy system and a sample, according to an embodiment of the invention.

FIG. 2 shows a system block diagram of a two-phase interferometric microscope system and a single emitter sample, according to an embodiment of the invention.

FIG. 3 shows a graph of the intensities of the radiated light detected at the two detectors of FIG. 2.

FIG. 4 shows a system block diagram of a two-phase interferometric microscope system and an extended sample, according to an embodiment of the invention.

FIG. 5 shows a graph of the intensities of the radiated light detected at the two detectors of FIG. 4.

FIG. 6 shows a system block diagram of a two-phase interferometric microscope system with representations of beam phases, according to an embodiment of the invention.

FIGS. 7A and 7B each show a graph displaying the phase of the beam intensities of the components for a given output beam detected at a respective detector of FIG. 6.

FIGS. 7C and 7D each show a graph displaying the intensities of the output beams associated with FIGS. 7A and 7B, respectively.

FIG. 8 shows a system block diagram of a three-phase interferometric microscope system, according to an embodiment of the invention.

FIG. 8A shows a system block diagram of a three-way beam splitter, according to an embodiment of the invention.

FIG. 8B shows a system block diagram of a eight-way beam splitter, according to yet another embodiment of the invention.

FIG. 9 shows a graph of the intensities of the radiated light detected at the three detectors of FIG. 8.

FIG. 10 shows a system block diagram of a three-way beam splitter, according to an embodiment of the invention.

FIGS. 10A, 10B and 10C show striped mirrors used as beam splitting devices in the three-way beam splitter of FIG. 10, according to an embodiment of the invention.

FIG. 11 shows a perspective view of the three-way beam splitter shown in FIG. 10.

FIG. 12 shows a system block diagram of a four-phase interferometric microscope system, according to an embodiment of the invention.

FIG. 13 shows a graph of the intensities of the radiated light detected at the four detectors of FIG. 12.

FIG. 14 shows a system block diagram of a four-phase interferometric microscope system having diffractive gratings, according to an embodiment of the invention.

FIG. 15 shows a system block diagram of a six-phase interferometric microscope system having diffractive gratings, according to an embodiment of the invention.

FIG. 16 shows a diagram of a single quantum emitter in two different potential locations.

FIG. 17 shows a system block diagram of a three-phase interferometric microscopy system having a single objective, according an embodiment of the invention.

FIG. 18 shows an example of a system block diagram of the beam splitter/annular reverser system shown in FIG. 17.

FIG. 19 shows a cross-sectional view of the annular expander shown in FIG. 18.

FIG. 20 shows an end view of the annular expander shown in FIG. 18.

FIG. 21 shows a perspective view of three out of the six prism components of the annular expander shown in FIG. 17.

FIG. 22 shows a graph of the intensities of the radiated light detected at the three detectors of FIG. 17.

FIG. 23 shows an end view of the annular expander, according to another embodiment of the invention.

FIG. 24 shows a cross-sectional view of the annular expander shown in FIG. 23 along line XX-XX.

FIG. 25 shows a cross-sectional view of the annular expander shown in FIG. 23 along line XXI-XXI.

FIG. 26 shows a flowchart of a calibration process for a multi-phase interferometric microscope system, according to an embodiment of the invention.

FIG. 27 shows a flowchart of a process for obtaining 3-dimensional position information of photo-activated samples within a multi-phase interferometric microscope system, according to an embodiment of the invention.

DETAILED DESCRIPTION

In broad terms, one or more embodiments described herein can produce 3-dimensional renderings with resolution better than known optical microscopy. In some embodiments, an interferometric microscope can measure a phase difference to determine position information in the third dimension (e.g., the z coordinate) while simultaneously measuring position information in the two other dimensions (e.g., the x and y coordinates). In some embodiments, an interferometric microscope can measure samples containing a collection of point-sources such as samples that internally radiate light from photo-activated fluorescent labels.

For example, 3-dimensional renders of specific labeled proteins in a sample can be produced. Contrast can be given by specific labeled proteins and still be acquired with a specialized optical microscope. For example, one or more embodiments can be used to locate a fluorescently labeled protein in a cell with full 3D location.

One known microscopy system is the Photo-Activated Localization Microscopy (“PALM”) system. This PALM system is described in the following references, which are incorporated herein by reference: U.S. Patent Application Ser. No. 60/683,337, entitled “Optical Microscopy with Phototransformable Optical Labels” and filed May 23, 2005; U.S. Patent Application Ser. No. 60/780,968, entitled “Imaging Intracellular Fluorescent Proteins at Near-Molecular Resolution” and filed Mar. 10, 2006; PCT Patent Application Serial No. PCT/US2006/019887, entitled “Optical Microscopy with Phototransformable Optical Labels,” filed May 23, 2006 and published on Nov. 11, 2005 as PCT Publication No. WO/2006/127682; Betzig, E. et al., “Imaging Intracellular Fluorescent Proteins at Nanometer Resolution,” Science, Vol. 313, Sep. 15, 2006; and Betzig, E. et al., “Supporting online Material for ‘Imaging Intracellular Fluorescent Proteins at Nanometer Resolution’,” Science Express, Aug. 10, 2006, [online] www.sciencemag.org/cgi/content/full/1127344/DC1 (collectively referred to herein as the “Betzig references”).

The PALM system localizes molecules predominantly in 2 dimensions, in the x,y image plane, to provide high-resolution images that exceed that of diffraction-limited microscopes. One known way to obtain information about the location of labeled molecules in the third dimension involves a serial process such as taking data on several sections, or by quantifying the point spread function variations of defocused molecules that depend on the distance to the focal plane.

As with 2D PALM system, switchable optical labels here can be the photo activatable fluorescent proteins such as PA-GFP. PA-GFP is a photo activatable fluorescent protein mutant of GFP where irradiation with 405 nm light can transform it from a non-fluorescent version into an anionic form that absorbs excitation radiation at 500 nm and emits fluorescent radiation at 520 nm. Other photo-activatable proteins are Dronpa, Eos, Kaede, Kikume, Kindling-FP, PA-CFP, many of which are useable in monomeric, dimeric and tetrameric forms. Such proteins can be genetically expressed and genetically attached to label a protein of interest. This enables protein specific images to be constructed. Various caged dye molecules can also be used but they typically require immuno-labeling or alternate attachment schemes. Because some of these proteins/labels have sufficiently different photophysical properties such as emission wavelength, excitation wavelength or bleaching rates, that it is possible to make images with two different switchable optical labels. In this case, the addition of dichroic mirrors in the beam paths can result in a new second or even third set of images for this 3D interferometric microscopy. Sequencing the acquisition to image and bleach one protein and then image the second is another way to image different labeled proteins.

Some embodiments described herein, however, use the phase angle information of the optical path length difference such that both two-dimensional intensity image, from which the x,y position of labeled molecules can be determined, and the phase dependent images, from which the z position of the labeled molecules can be determined, can be simultaneously acquired and can be reconstructed. In other words, two or more emission beams from a sample are further split into two or more beams and pairwise interfered with each other (based on their differing phases) and the resulting interfered beams form images that are recorded simultaneously to reconstruct the phase of the path length difference (i.e., the z position) of the localized sources of a sample. In addition, in some embodiments, the interference amplitude can also be extracted and used as a supplemental measure of vertical position.

Such embodiments are advantageous in many respects. For example, as mentioned above, in some embodiments described herein, the phase difference can be measured to determine position information in the third dimension (e.g., the z coordinate) while simultaneously measuring position information in the two other dimensions (e.g., the x and y coordinates). This allows for a measurement of all three dimensions that can be more accurate than known systems where these measurements occur at different times. Such a measurement technique described herein can be particularly advantageous, for example, for samples that are not static such as photo-activated samples that have a time-varying brightness.

As also mentioned above, in some embodiments described herein, the component of the sample being measured is separable and identifiable as a sub-diffractive sized point source. This allows for a measurement that can be well defined and more accurate than known systems where the measured sample component is an extended source with a distribution of positions and various portions of the source can interfere with other portions of the source. Such a measurement technique described herein can be particularly advantageous for performing measurements with an accuracy below the diffraction limit.

FIG. 1 shows a system block diagram of a two-phase interferometric microscopy system and a sample, according to an embodiment of the invention. As shown in FIG. 1, the microscopy system 100 includes optical components 130 and 140, a beam splitter 150, optical components 160 and 170, and detectors 165 and 175.

The optical components 130 and 140 can be any type of optical component(s) that define an optical path 132 and 142, respectively, from the sample S to the beam splitter 150. For example, optical components 130 and 140 each can include a set of mirrors, lenses or objectives. The optical components 130 and 140, for example, can include lenses or objectives that collimate optical energy emitted from the point source of the sample. Although FIG. 1 shows two diametrically opposite paths 132 and 142 into two optical components 130 and 140, in some embodiments discussed below, a single optical component such as optical component 140 can encompass two non-opposing paths 132 and 142.

Optical components 130 and/or 140 can include additional components and/or features that allow the sample S to receive one or more sources of activation. For example, in some embodiments, optical energy from optical sources (not shown in FIG. 1) can be used to activate the fluorescent response of the sample S to excitation photons. In such embodiments, for example, the optical components 130 and/or 140 can include an aperture through which optical energy from the optical sources can irradiate the sample S while minimizing stray radiation from the excitation laser traveling back towards the detectors 165 and/or 175. Laser line filters can also be in the beam path to the detectors 165 and/or 175 so that predominantly fluorescent light reaches the detectors. A first laser source can radiate the sample S, for example, with an activation beam and a second laser source can radiate the sample S, for example, with an excitation beam. In other embodiments, an optical source other than a laser, such as a light-emitting diode, incandescent lamp, can be used as the source for activation and/or excitation. As described in further detail in the Betzig references incorporated by reference above, the activation beam is configured to activate sufficiently sparse subsets of individually resolvable photo-activatable labels in the sample, and the excitation beam is configured to excite the activated portions of the sample S to emit fluorescence-based photons.

In alternative embodiments, the fluorescent response of the sample can be activated with one or more chemicals. In such alternative embodiments, the sample can be exposed to a first chemical to activate the fluorescent sparse subsets of molecules of the sample and subsequently to a second chemical to deactivate the fluorescence of the sample if photo-bleaching is not used.

The beam splitter 150 can be any type of appropriate optical component that receives and combines the optical energy along optical paths 132 and 142, and sends the combined and mixed optical energy along optical paths 152 and 154 (each also referred to herein as “a leg”) to optical components 160 and 170. For example, beam splitter 150 can be a 50:50 beam splitter where the combined optical energy is split such that substantially 50 percent of the combined optical energy is sent on optical path 152 and substantially 50 percent of the combined optical energy is sent on optical path 154. As optical energy from optical paths 132 and 142 is interfered and mixed by the beam splitter 150, the resulting optical energy along optical path 152 is based on a relative phase difference between the optical energies from optical paths 132 and 142. Similarly, as optical energy along optical paths 132 and 142 is interfered and mixed in the beam splitter 150, the resulting optical energy along optical path 154 is based on a relative phase difference between the optical energies from optical paths 132 and 142, which is different from the relative phase difference associated with the optical energy along optical path 152. In this embodiment, the resulting optical energy along optical path 152 and optical energy along optical path 154 differ in phase by approximately 180 degrees. Alternatively, beam splitter 150 can divide the beam into any desired proportions such as 66:33, 70:30, etc. In yet another embodiments, beam splitter 150 can be a combination of multiple beam splitters. Such beam splitters can be constructed, for example, from a block of two glass segments and a thin film coating disposed between the two glass segments. Such beam splitters could produce three or more outputs, i.e., more than the two of 152 and 154. This enables three or more output beams to be expressed with three or more interference phase angles (e.g., 0 degrees, 120 degrees and 240 degrees), rather than 0 and 180 degrees in the two beam cases. In yet other embodiments, in addition to or alternative to a beam splitter having a thin film coating, diffractive gratings can be used as beam splitters and can divide a beam into two or more beams. Gratings can also take multiple input beams and mix and interfere them into various output beams. In some embodiments, the diffraction gratings can be used in conjunction with a sample having point sources emitting at the substantially same wavelength such as a label having a switchable scattering center (e.g., a photoswitchable label with an activatable scattering characteristic). Such a sample having point sources emitting at substantially the same wavelength advantageously allows the diffraction gratings to deflect input beams at a relatively narrow deflection angle, which is a function of the wavelength of the input beam.

In another embodiment, one or more of the beam splitters can include striped mirror portions instead of diffraction gratings or thin-film coatings. In one embodiment, the striped mirror portions have a stripe period that is sufficiently small such that the diffractive order does not refocus into the image plane. For a wavelength λ, focal length f, and image size d, this means a period less than fλ/d. Other patterns could be used such as square or hexagonal polka dots. To prevent ghosting in the image, however, the same condition of the polka dot period (i.e., less than fλ/d) must be maintained.

Similar to optical components 130 and 140, optical components 160 and 170 can be any type of optical component(s) that defines an optical path 152 and 154, respectively, from the beam splitter 150 to the detectors 165 and 175, respectively. For example, optical components 160 and 170 can be lenses that de-collimate optical energy from optical paths 152 and 154, respectively, so that the optical energy is focused at detectors 165 and 175, respectively. In some embodiments, one or more of the optical components such as optical components 130, 140, 160 and 170, and/or beam splitter 150 can include adaptive optics.

In use, a statistically sparse subset of separately-imaged labels in the sample S is activated before a measurement is made. As mentioned above, the statistically sparse subset of the labels in the sample S can be activated by an appropriate technique such as by photo activation or chemical activation. Once activated, the sample S can be excited by an excitation source such as a laser. This results in photons being emitted from the sample S along optical paths 132 and 142 through optical components 130 and 140 to beam splitter 150. The optical energy from optical paths 132 and 142 are combined at beam splitter 150 and sent along paths 152 and 154 through optical components 160 and 170 to detectors 165 and 175, respectively. A phase difference between the beams 132 and 142 causes interference and modulates the intensity of the detected signals at detectors 165 and 175, which can used to determine a position in the z coordinate, as described below in further detail. In addition, the detected signals at detectors 165 and 175 also provide position information in the x, y coordinates. After the detected signals are collected, the activated sparse subsets of the labels in the sample S can be bleached by repeated excitation so that no further optical signals can be detected from those labels, and further sparse subsets of labels can then be activated and excited for additional measurements, and the process repeated.

FIG. 2 shows a system block diagram of a two-phase interferometric microscopy system and a sample with a single active emitter, according to an embodiment of the invention. As shown in FIG. 2, a single quantum emitter S is located a distance, δ, from a center location 220. The single quantum emitter S can be, for example, a fluorescent molecule, a quantum dot or a sub-wavelength scattering center. A radiated photon can escape from the quantum emitter S both upwards into lens 230 and downwards into the lens 240, where two collimated beams 232 and 242, respectively, are formed. The collimated beams 232 and 242 reflected by mirrors 234 and 244, respectively, and are received at beam splitter 250.

Two beam paths are defined: a first beam path (also referred to as the top path) from the quantum emitter S to mirror 234 and from mirror 234 to beam splitter 250; and a second beam path (also referred to as the bottom path) from the quantum emitter S to mirror 244 and from mirror 244 to beam splitter 250. The mirrors 234 and 244 and beam splitter 250 are positioned relative to the quantum emitter S such that the lengths of the beam paths are matched to about a wavelength (L₁+L₂) top=(L₁+L₂) _(bottom), with the exception of a positional displacement δ of the emitter and a 90 degree phase shift on the top beam path, where L₁ is the path length between the center location 220 and the mirrors 234 or 244, and L₂ is the path length between the mirrors 234 or 244 and the beam splitter 250. In this case the electric field of the optical beam in the top beam path is: E ₁ =A*exp(−ikδ+ikL ₁ +ikL ₂ +iπ/2),

where k=2π/λ.

The electric field of the optical beam in the bottom beam path is: E ₂ =A*exp(ikδ+ikL ₁ +ikL ₂).

When received at the beam splitter 250, the optical beams are combined and transmitted to the two detectors 265 and 275 through lenses 260 and 270, respectively. The electric field from one source of the optical beam going towards each detector 265 and 270 is: E ₁₂=½^(0.5)*(E ₁ ±E ₂)=½^(0.5)*(A*exp(−ikδ+ikL ₁ +ikL ₂ +iπ/2)±A*exp(ikδ+ikL ₁ +ikL ₂)).

In this embodiment, the detectors 265 and 275 each is a charged coupled device (CCD) detector or (for improved sensitivity to single fluorescent molecules) an electron-multiplying charged coupled detector (EMCCD). As shown in FIG. 3, the intensity 365 of the radiated light detected at detector 265 and the intensity 375 of the radiated light detected at detector 275 vary with a sine pattern with phase 4πδ/λ. More specifically, the intensity 365 of the radiated light detected at detector 265 is: I ₁=½*(E ₁ +E ₂)² =A ²*(1+sin(2kδ)).

The intensity 275 of the radiated light detected at detector 275 is: I ₂=½*(E ₁ −E ₂)² =A ²*(1+sin(2kδ+π)).

The displacement can then be calculated by combining and inverting the two intensity equations: (I ₁ −I ₂)/(I ₁ +I ₂)=sin(4πδ/λ) or δ=λ arcsin((I ₁ −I ₂)/(I ₁ +I ₂))/4π

This can be generalized to a microscopy system such as that shown in FIG. 4. More specifically, FIG. 4 shows a system block diagram of a two-phase interferometric microscopy system and an extended sample, according to an embodiment of the invention. As shown in FIG. 4, the microscopy system includes microscope objectives 430 and 440 and lenses 460 and 470 such that an image plane is defined at detectors 465 and 475. Any point source at sample S will form a resolution limited spot on each image plane defined at detectors 465 and 475.

The intensity of the spatially integrated spot can then be used to calculate the z-coordinate position of the point source of the sample using an intensity equation described above. By using a CCD, EMCCD, or other large area detector, multiple image spots can be detected simultaneously and the z-displacement δ of each associated portion of the sample 310 can be established at substantially the same time. The phase of the detected optical signals can be determined to a fraction of a radian, so that positional information can be determined to better than a fraction (for example, <20%) of a wavelength. For example, for a wavelength of λ=600 nm, the positional information can be determined to less than 1 radian, which corresponds to δ=50 nm. Such an interferometric system can be configured such that the x,y locations of the source emitters can also be measured to the nanometer scale resolution.

The thickness of the sample can affect the extent to which a unique z-coordinate position can be determined. For example, in the two-phase embodiments shown in FIG. 2 or 4, if the thickness of the sample 210 is less than λ/4, then a unique z-coordinate position can be determined. The sensitivity to the displacement, δ, for these embodiments is non-uniform. In other words, if the phase is such that intensity is determined close to the minimum or maximum of the sine wave, then the slope approaches zero and therefore the sensitivity approaches to zero as shown in FIGS. 3 and 5. If, however, the phase is such that the intensity is determined away from the minimum or maximum of the sine wave, then the sensitivity to the displacement, δ, for these embodiments is measurable and nonzero.

FIG. 6 shows a system block diagram of a two-phase interferometric microscope system with representations of beam phases, according to an embodiment of the invention. Beam 432 has a phase different from the phase of beam 442 and output beam 452 includes components 452′ and 452″ which differ in phase. Although beam components 452′ and 452″ are shown as being separated in FIG. 6, this is for illustration purposes and it should be understood that the beam components are coincidental and interfering with each other. Beam component 452′ originates from beam 432 and beam component 452″ originates from beam 442. Similarly, output beam 454 includes components 454′ and 454″ which differ in phase. Although beam components 454′ and 454″ are show as being separated in FIG. 6, this is for illustration purposes and it should be understood that the beam components are coincidental and interfering with each other. Beam component 454′ originates from beam 432 and beam component 454″ originates from beam 442. In this embodiment, the phase difference between beam components 452′ and 452″ determines an intensity that is detected at detector CCD φ₁ and the phase difference between beam components 454′ and 454″ determines an intensity that is detected at detector CCD φ₂.

FIGS. 7A and 7B each show a graph displaying the phase of the beam intensities of the components for a given output beam detected at a respective detector of FIG. 6. Output beams 452 and 454 in FIG. 6 include components of input beams 432 and 442. As shown in FIG. 7A, 452′ and 452″ are components of output beam 452 and are separated in phase by a phase difference 480. As shown in FIG. 7B, 454′ and 454″ are components of output beam 454 and are separated in phase by a phase difference 490. The difference in these phase differences will vary as a function of the position of the object. The squared sum of the beams represented in FIG. 7A produce a beam with the intensity shown in FIG. 7C, and likewise the squared sum of the beams in represented FIG. 7B produce a beam with the intensity shown in FIG. 7D. The relative intensity of these two beams, I₁ and I₂, determine the interference phase and thereby the vertical displacement δ.

In alternative embodiments, the interferometric microscopy system can have more than two-phases. For example, FIG. 8 shows a system block diagram of a three-phase interferometric microscopy system, according to an embodiment of the invention. FIG. 9 shows a graph of the intensities of the radiated light from a point source when properly interfered and detected at the three detectors of FIG. 8.

FIGS. 10 and 11 show an embodiment of a 3-way beam splitter. The two emission beams 1010 and 1020 enter a 50-50 beam splitting surface 1050. The two resulting interfered beams proceed to an upper 66:33 beam splitting surface 1040 and a lower mirror surface 1070. Most of the upper beam, 66%, exits into output path 1030 and 33% goes to path 1028 and interferes with the beam 1024 that has been reflected from the mirrored surface 1070. This final interference takes place on a 50:50 beam splitting surface 1052 resulting in the second and third output beams 1032 and 1034. To have a balanced 3-way interference (i.e., similar amplitudes for output beams 1030, 1032 and 1034) between the 3 outputs, the phase of the interfering beams 1028 and 1024 can be tuned by changing the position of the mirror 1070 and thereby the gap 1060. The gap 1060 can be filled with index matching oil to minimize residual reflectivity on the bottom of the glass prism.

FIGS. 10A, 10B and 10C show striped mirrors used as beam splitting devices in the three-way beam splitter of FIG. 10, according to an embodiment of the invention. Each striped mirror shown in FIGS. 10A, 10B and 10C is rotated with respect to the other striped mirrors shown in FIGS. 10A, 10B and 10C to reduce undesirable interference in the beam paths.

FIGS. 10 and 11 further represent propagation of wavefronts associated with object and image planes through an embodiment of a 3-way beam splitter. Emission beams 1010 and 1020 having wavefronts 1012 and 1022 enter beam splitting surface 1050. Wavefronts 1012 and 1022 are associated with the object plane of the sample and are interfered with each other at beam splitting surface 1050. The resulting wavefronts interfere with each other at beam splitting surfaces 1040 and 1052, as they propagate through the 3-way beam splitter. Beam splitting surfaces 1040, 1050 and 1052 have sufficient surface area to accommodate distortion or spreading of the wavefronts during interference at beam splitting surfaces 1040, 1050 and 1052. Output beams 1030, 1032 and 1034 emerge from the 3-way beam splitter having combined the beams 1010 and 1020 with a relative phase shift of 1036, 1038 and 1039. The results in different intensities of beams 1030, 1032 and 103, which are intact, interfered versions of wavefronts 1012 and 1022 and propagate to the image planes of the detectors.

FIG. 12 shows a system block diagram of a four-phase interferometric microscopy system, according to an embodiment of the invention. FIG. 13 shows a graph of the intensities of the radiated light detected at the four detectors of FIG. 12. Although FIGS. 8 through 12 show embodiments having three- or four-phase interferometric microscopy systems, systems having more than four phases are possible such as five-phase interferometric microscopy systems, six-phases interferometric microscopy systems, etc.

While a two-phase interferometric microscopy system can use, for example, a beam splitter such as 50:50 beam splitter 250 as shown FIG. 2, the beam splitter for microscopy systems having more than two phases (e.g., as shown in FIG. 8) can be embodied in a number of different ways. For example, as shown in FIG. 8, the beam splitter can include multiple components: two 50:50 beam splitters, a 66:33 beam splitter and a mirror. These beam splitter components can be positioned manually, for example, on an optical bench. Alternatively, these beam splitter components can be manufactured and disposed between a pair of optical flats; the position of such beam splitter components can be adjusted, for example, by piezoelectric transducers. Such adjustments can insure that the phase shifts of the three legs are spaced by 2π/3 or for a N leg embodiment by 2π/N and also parallelness of the splitting planes is established. An alternative embodiment of a 3-way beam splitter is shown in FIG. 8A.

FIG. 8B shows a system block diagram of a eight-way beam splitter, according to yet another embodiment of the invention. As shown in FIG. 8B, the beam splitter includes multiple beam splitter components each of which is represented by a dashed line. Although FIG. 8B explicitly shows a beam splitter with 16 beam splitter components that produces 8 output beams, it should be understood that the basic structure can be followed for a greater or lesser number of beam splitter components to produce N number of output beams. These beam splitter components can have a reflectivity that ranges between 0 to 100%. Transmission can be nominal. When the surfaces of the beam splitter components are substantially parallel to each other, substantially the whole extent of a beam (or beam portion) (having a diameter, D) can interfere coherently with reflected beam portions of or other beam portions split from the original input beams A and B. The surfaces of the beam splitter components can be substantially parallel, for example, with each being within an angle of <λ/D. Also the beam splitter components are equally spaced to within a coherence length λ²/(nΔλ) of the input beams, where λ is wavelength, Δλ is a range of wavelengths, and n is the index of refraction of the medium. This spacing insures that the interference is substantially coherent across the beam wavelength range dl of interest. Input beams A and B enter the beam splitter with a direction vector so that A-B is a vector normal to the planes of the beam splitter. Both input beams are aimed to the same point on the first beam splitter component. The parallel planar structure of the beam splitter allows a range of beams angles to self interfere; this range of beam angles can correspond to the range of angles associated with an image in a collimated portion of a beam before it is focused onto an image plane of, for example say a CCD detector (not shown in FIG. 8B). In some embodiments, both the reflectivity of each beam splitter component and the fine spacing between beam splitter components (less than the coherence length but enough to control interference angle) can be varied to ensure a substantially equal splitting into the N output beams (FIG. 8B shows N=8) for the case where only one input beam (e.g., input beam A or input beam B) is applied. Furthermore, in some embodiments, the interference angle between the input beam A and the input beam B can be set to different intervals of 360 degrees/N. So, in the example shown in FIG. 8B where N=8, the intervals of 45 degrees, 90 degrees, 135.degrees can be established.

As mentioned above, in some embodiments, one or more beam splitters can be based on diffraction gratings instead of beam splitters having thin-film coatings. For example, FIG. 14 shows a system block diagram of a four-phase interferometric microscope system having diffractive gratings, according to an embodiment of the invention. In this embodiment, the beam splitter 1450 includes diffraction gratings 1451 through 1454, mirrors 1462 and 1464, and phase shifters 1472 and 1474. Diffraction gratings 1451 and 1452 receive optical energy from mirrors 1434 and 1444, respectively. Diffraction grating 1451 divides the received optical energy and sends the divided portions to mirror 1462 and phase shifter 1472. Diffraction grating 1452 divides the received optical energy and sends the divided portions it to mirror 1464 and phase shifter 1472. Diffraction grating 1453 receives and mixes the optical energy from mirror 1462 and phase shifter 1474, divides the mixed optical energy, and sends each of the divided portions to a different detector. Similarly, diffraction grating 1454 receives and mixes the optical energy from mirror 1464 and phase shifter 1472, divides the mixed optical energy, and sends each of the divided portions to a different detector. FIG. 15 shows a system block diagram of a six-phase interferometric microscope system having diffractive gratings, according to another embodiment.

In some embodiments, spectral filters (e.g., bandwidth filters, low-pass filters or high-pass filters) can be disposed in the optical paths to reduce or prevent particular optical energies from contributing to the optical energy emitted from the switchable optical source. For example, optical energy used during photo-activation and/or photo-excitation can be removed from the optical paths to decrease stray light and noise within the optical paths and propagating towards the detectors.

In some embodiments, one or more the sample, objectives, mirrors, optical components, beam splitters, lenses, spectral filters and other elements along the optical paths can be coupled individually to multi-axis positioners including, for example, piezoelectric actuators. Such multi-axis positioners can allow for automated or non-automated calibration and positioning. Additionally, such multi-axis positioners can improve precision and accuracy in positioning and calibration over manual methods.

During setup and as a guide to positioning and alignment of the interferometer system, it can be useful to have a shutter incorporated into the two input beam paths. This allows one to observe the image from each path independently, and make sure that each path is optimally focused. Furthermore, because the images of the two portions can be observed, it is possible to shift one image with respect to the other so that they overlayed for optimal interference when the shutters are open for both beams.

In alternative embodiments having more than two-phases, the sensitivity can be phase invariant. In addition, the unique range of operation can be extended to λ/2. Larger displacements will result in a positional ambiguity modulo λ/2. This can be lifted with supplemental techniques, such as, for example, selective illumination (excitation or activation) in both space and time domain, or sample sectioning. In such embodiments, for example, having three-phase or four-phase interferometry, the position of the measured portion of the sample is determined by the intensities: I _(ij) =G _(j)*(A _(Ui) ² +A _(Li) ²+2*β_(j) *A _(Ui) *A _(Li)*cos((φ_(i)+φ_(j))+Ofs _(j))

where j=0, 1, 2, . . . indicates the leg of multiphase interferometer;

I_(j)=intensity of the j leg of point source, with offset no correction;

β_(j)=interference correction factor of j^(th) leg;

G_(j)=gain of the j^(th) leg;

Ofs_(j)=offset of the j^(th) leg;

A_(Ui)=Amplitude of i^(th) point source upper objective;

A_(Li)=Amplitude of i^(th) point source lower objective;

φ_(j)=j^(th) leg phase shift; and

φ_(i)=point source (i.e., fluorescent molecule) phase.

Dark calibration A_(Ui)=A_(Li)=0 allows one to determine the dark signal from the j^(th) detector I _(oj) =G _(j) *Ofs _(j) I _(ij) =G _(j)*(A _(Ui) ² +A _(Li) ²+2*β_(j) *A _(Ui) *A _(Li)*cos(φ_(i)+φ_(j))+Ofs _(j))Total amplitude: A _(Ti) ² =A _(Ui) ² +A _(Li) ²

Rescaled interfering fraction: X _(i)=2*A _(Ui) *A _(Li) /A _(Ti) ²

So, with these variables the intensities are: I _(ij) =I _(oj) +G _(j) *A _(Ti) ²(1+β_(j) *X _(i)*cos(φ_(i)+φ_(j)))

A calibration where a fiducial source is moved by δ along z can be used to determine the remaining apparatus parameters G_(j), β_(j) and φ_(j). This will result in an intensity signal that oscillates in each leg similar to FIG. 7 for the three phase case as the phase φ_(f)(z)˜4*π/(nλ)δ(z) winds. Here, n is the index of media with the shifted path length (e.g., immersion oil).

A fit to set of such calibration data will then determine the following equation for each leg j on a reference fiducial with interference fraction of the emitting fiducial X_(f) and its intensity A_(Ti) ²=I_(f) where the apparatus parameters G_(j), β_(j) and φ_(j) are now determined. Such calibration can be useful in the phase-shift values are no exactly 120 degrees (3 legs) or 2π/N (N legs). I _(fj)(z)=I _(oj) +G _(j) *I _(f)(1+β_(j) *X _(f)*cos(φ_(i)(z)+φ_(j))).

G_(j), β_(j) can be defined by normalizing to fiducial units: I _(f)=1 and X_(f)=1 and α_(i) =A _(Ti) ² /I _(f) I _(ij) =I _(oj) +G _(j) *A _(Ti) ²(1+β_(j) *X _(i)*cos(φ_(i)+φ_(j))) D _(ij)=(I _(ij) −I _(oj))/(G _(j) *I _(f)) rescaled intensity data D _(ij)=α_(i)*(1+β_(j) *X _(i)*cos(φ_(i)+φ_(j)))

FIG. 26 shows a flowchart of a calibration process 2600 for a multi-phase interferometric microscope system using an automated calibration method, according to an embodiment of the invention. At 2610, an optical emitter including sub-diffractive fiducials is imaged at various positions in the z plane to produce an image frame sequence. This process at 2610 includes sub-processes 2611 and 2612, which are repeatedly performed for various positions in the z plane to produce the image frame sequence. At 2611, an image of the optical emitter is captured at each detector in the interferometric microscope system. The image frames captured at each detector contain phase-related characteristics (e.g., φ_(j) β_(j) of the j^(th) detectors) that are phase-shifted relative to the information in the image frames captured at the other detectors. The z position of the optical emitter is incremented at 2612 and 2611 is repeated N times to produce an image frame sequence. Each image frame in an image frame sequence is based on optical energy intensities in the x,y plane relating to the emission distribution of the optical emitter.

At 2621, the emission distributions of the optical emitter are fit to a predetermined distribution to localize the optical emitter emission distribution in the x,y plane. For example, certain optical emitters have emission distributions that can be accurately fit to Gaussian distributions. Such Gaussian distributions can represent, for example, a likelihood a point spread function of the images of optical emitters. At 2622, the image frames are corrected for drift using the fiducials included in the image frames. The drift-corrected frames are then transformed at 2623 to produce transformed frame that align the fiducials in the image frames in the x,y plane.

At 2624, the image frames captured at each of the image detectors for each position in the z plane are summed into summed image frames. The summed image frames are a composite of the phase-shifted image frames captured by each detector. At 2625, the summed image frames are fit to a predetermined distribution to determine the center of the optical emitter in the x,y plane.

The transformed frames are fit at 2626 to a predetermined distribution, such as a Gaussian distribution, and the amplitude of the optical emitter emission distribution in each frame is measured at the center of the optical emitter in the x,y plane determined in 2625. The center of the optical emitter can be based on, for example, the centroid of a Gaussian distribution (e.g., a three-dimensional distribution) representing a likelihood of the location of an optical emitter. At 2627, the amplitudes of the optical emitter emission distribution are normalized to the sum determined at 2624. The normalized amplitudes determined at 2627 for each image and the associated z position is recorded at 2628. For example, the normalized amplitudes and the associated z positions for each frame can be stored within a look-up table or parameterized by the previous equations to determine φ_(j) and β_(j).

The following equations set forth an explanation by which the displacement, δ, of source (i.e., the relative z-coordinate position of the sample) can be determined from the intensity data collected by the interferometry microscopy system.

Once the calibration is complete, the data can be scaled and then the values β_(j) and φ_(j) are also known for each output beam of the interferometer. The j equations: D _(ij)=α_(i)*(1+β_(j) *X _(i)*cos(φ_(i)(z)+φ_(j))) can be used to determine the position of an emitter along the axial direction. Thus three or more data values D_(i0), D_(i1), D_(i2), . . . for the intensity of the i^(th) molecule result in three or more equations to determine the three unknown values: α_(i), X_(i), and (of most interest) φ_(i)(z). If there are only three phases, then the three equations will yield a unique solution for these three unknowns. This and the value of φ_(i) in particular can be solved for example by Newton's Method.

If there are more than three phases, then four or more equations can determine the three unknown values and the system of equations is over constrained. In that case, the best values for the unknowns can be determined using a chi-squared minimization method. By minimizing the error of the equations: χ_(i)=Σ_(j) [D _(ij)−α_(i)(1+β_(j) cos(φ_(i)+φ_(j)))]²

In the approximation of a well balanced interferometer system (meaning phase shift increment between cameras and equalized gain and equalized interference efficiency) we can obtain a best value for φ_(i) given by the equation: tan(φ_(i))≅−Σ_(j) D _(ij) sin(φ_(j))/Σ_(j) D _(ij) cos(φ_(j))

Once the phase φ_(i) of the i^(th) emitter is known then the vertical position δ_(i) along the z axis can be deduced to a first approximation by the equation. Note n is the index of refraction of the sample environment. δ_(i)≅λ*φ_(i)/(4πn)=λ/(4πn)*atan(Σ_(j) D _(ij) sin(φ_(j))/Σ_(j) D _(ij) cos(φ_(j)))

This last equation solving for displacement, δ, describes the phase shift from the intensity data. In other words, this last equation shows how the displacement, δ, can be calculated based on the intensity data detected at the detector for each leg of the interferometry microscope system. For example, for a 3-phase interferometry microscope system (φ_(j))≅−2π/3, 0, 2π/3. For a 4-phase interferometry microscope, (φ_(j))≅0, π/2, π, 3π/2.

Note that the above equation for calculating displacement, δ, is based on various assumptions. In alternative embodiments, other assumptions can be used resulting in one or more alternatives for calculating the displacement, δ. For example, higher order corrections can be calculated from the real part and taking out the approximations that the N phases of each interferometer leg are equally spaced or that the sensitivity of interferometer leg is equal. Alternative embodiments can include further refinements that provide a means to calibrate and correct for non-balanced intensity sensitivities and phase mixing angles.

In summary, the N spot intensities of a N phase interferometer becomes the data D_(j) that can be converted to the vertical position displacement δ.

Because φ_(i) is a cyclic solution, for every 2π the uniquely determinable range extends λ/2n. For samples thicker than this, an ambiguity exists in determining δ by +/−N λ/2n for an integer N. Other attributes, such as measured point spread function and interference efficiency, can be used to find the most likely N and chose the best of a multi-valued solution. In effect the extent of defocus is an alternate, less precise but adequate method to estimate to within λ/2n where the emitter is located. This combined with the interferometricly-deduced position can give a unique vertical position of the emitter. One further consideration in estimating the vertical position is to account for the Gouy phase, an extra term that advances the phase by an extra π at the focus over about a Rayleigh length distance z_(R). In this case δ is given by the solving the following equation: φ_(i)≅4πδ/(nλ)+atan(δ/z _(R)).

FIG. 27 shows a flow chart of a process for obtaining 3-dimensional position information of photo-activated samples within a multi-phase interferometric microscope system, according to an embodiment of the invention. At 2710, various switchable optical sources including sub-diffractive fiducials are imaged to produce an image frame sequence. This process at 2710 includes sub-processes 2711 and 2712, which are repeatedly performed for various positions in the z plane to produce the image frame sequence. At 2711, a switchable optical source is activated and excited. At 2712, an image frame is captured at each detector in the interferometric microscope system and 2711 is repeated N times to produce an image frame sequence. The image frames captured at each detector contain phase-related characteristics that are phase-shifted relative to the information in the image frames captured at the other detectors. The image frames are based on optical energy intensities in the x,y plane relating to the emission distributions of the switchable optical sources.

At 2721, the emission distributions of the switchable optical sources are fit to a predetermined distribution to localize the emission distributions of the switchable optical sources in the x,y plane. For example, certain switchable optical sources have emission distributions that can be accurately fit to Gaussian distributions. Such Gaussian distributions can represent, for example, a point spread function of the images of optical emitters. At 2722, the image frames are corrected for drift using the fiducials included in the image frames. The drift-corrected frames are then transformed at 2723 to produce transformed frames that align the fiducials in the image frames of the image frame sequence in the x,y plane.

At 2724, the image frames captured at each of the image detectors for each position in the z plane are summed into summed image frames. The summed image frames are a composite of the phase-shifted image frames captured by each detector. At 2725, the summed image frames are fit to a predetermined distribution to determine the centers of the switchable optical sources in the x,y plane.

The transformed frames are fit at 2726 to a predetermined distribution, such as a Gaussian distribution, and the amplitude of the emission distribution of the switchable optical source in each frame is measured at the centers of the switchable optical sources in the x,y plane determined in 2725. The center of the optical emitter can be based on, for example, the centroid of a Gaussian distribution representing a likelihood of the location of an optical emitter. At 2727, the amplitudes of the emission dispersions of the switchable optical source are normalized to the sum determined at 2724. At 2728, the normalized amplitudes determined at 2727 are used to obtain the z plane coordinate of the switchable optical source by correlation with calibration data recorded at previous step 2528. Such calibration data can be stored, for example, in a lookup table. The z plane coordinated can be calculated based on the normalized amplitudes, for example, by determining the minimum chi-square of: [CA₁(z)−PA₁]²+[CA₂(z)−PA₂]²+[CA₃(z)−PA₃]²+ . . . +[CA_(N)(z)−PA_(N)]²

-   -   where N=the number of detectors in the system;     -   CA_(i)(z)=the normalized peak optical energy intensity amplitude         at detector _(i) determined during calibration for a particular         sample position in the z plane;     -   PA_(i)=normalized peak optical energy intensity amplitude at         detector _(i) determined from the current PALM image frame.

In another embodiment, an interferometric microscopy system can use different portions of the emitted beam that still generate a path length difference that varies with a third dimension z. Such an alternative embodiment is discussed in connection with FIGS. 16-22. More specifically, FIG. 16 shows a diagram of a single quantum emitter in two different potential locations; FIGS. 17 and 18 show an example of a system block diagram of a three-phase interferometric microscopy system having a single objective; FIGS. 19-21 show an example of an annular expander of FIG. 18; FIG. 22 shows a graph of the intensities of the radiated light detected at the three detectors of FIG. 17. These figures are discussed below.

FIG. 16 shows an example of a single quantum emitter, in two possible locations. Such a quantum emitter can be, for example, a fluorescent molecule, a quantum dot or a sub-wavelength scattering center. A radiated photon (wavelength λ) will propagate in all directions including a downward direction away from the objective shown in FIG. 17 and an upward direction into the objective shown in FIG. 17. The portion of the energy of the radiated photon propagating into the objective can be considered as two collimated beams of roughly equal strength: an inner cylindrical beam labeled “B” and an outer annular beam labeled “C.” The size of inner beam portion B and outer beam portion C can be selected such that their intensities are substantially equal; such substantially equal intensities correspond to different beam portion areas due to the non-linear intensity profile of the beam (e.g., the peak of the beam energy is in the center of the beam with the intensity reducing towards the edges of the beam). For example, the inner beam portion labeled “B” can have a radius equal to approximately one third of the total beam radius; the outer beam portion labeled “C” can have an inner radius corresponding to the radius of the inner beam portion “B” and an outer radius corresponding to the radius of the total beam.

When the sample is exactly in focus, these beam portions are matched and in phase (as shown in the left side of FIG. 16). When, however, there is a positional displacement Δz of the emitter (as shown in the right side of FIG. 16), the phase beam portions B and C are displaced by Δ_(inner) and Δ_(outer), respectively, where Δ_(inner) equals Δz and Δ_(outer) equals Δz*cos(θ). This displacement of the phase beam portions B and C is due to the fact that the outer beam portion C will undergo less of an optical path change by a factor of cos(θ), where θ is the initial cone angle of rays from the emitter. Consequently, a differential phase shift between the two beams remains and equals: δ=2π(Δ_(inner)−Δ_(outer))/λ=Δz*(1−cos(θ))*2πn/(λ).

The z range over which the phase increases through 360 degrees and repeats is: λ/(n(1−cos(θ)). This z range for this single-objective embodiment is longer than the z range of the opposing-objective embodiments discussed above in connection with FIGS. 1-11. The index of refraction, n, of the optical environment around the emitter rescales the phase shift. An embodiment of a single-objective system is discussed in more detail below in connection with FIGS. 17-20.

FIG. 17 shows an example of a system block diagram of a three-phase interferometric microscopy system having a single objective. FIG. 18 shows an example of a system block diagram of the beam splitter/annular reverser system shown in FIG. 17. As shown in FIG. 17, the interferometric microscopy system includes an objective 1730; a mirror 1744; a beam splitter/annular reverser system 1740; mirror 1756; beam splitters 1754, 1755 and 1758; lens 1760, 1770 and 1780; and detectors 1765, 1775 and 1785. As shown in FIG. 18, the beam splitter/annular reverser system 1740 includes a beam splitter 1741, mirrors 1742-1747, annular expander 1748 and annular reducer 1749.

The annular expander 1748 can be any type of device that receives light having a given radius and expands the radius of the light. For example, the annular expander 1748 can receive the inner cylindrical beam of light from the mirror 1742 (similar in cross-sectional shape to the inner beam “B” discussed in reference to FIG. 17) and can produce light having an annular shape with a larger outer diameter (e.g., similar in cross-sectional shape to the outer beam “C” discussed in reference to FIG. 17). Similarly, annular reducer 1749 can be any type of device that receives light having a given radius and reduces the radius of light. For example, the annular reducer 1749 can receive the outer annular beam of light from the mirror 1742 (similar in cross-sectional shape to the outer beam “C” discussed in reference to FIG. 17) and can produce light having an inner cylindrical shape (or smaller annular shape) with a smaller inner diameter (e.g., similar in cross-sectional shape to the inner beam “B” discussed in reference to FIG. 17).

FIGS. 19-21 show an example of an annular expander shown in FIG. 18. More specifically, FIG. 19 shows a cross-sectional view of the annular expander shown in FIG. 18; FIG. 20 shows an end view of the annular expander shown in FIG. 18; and FIG. 21 shows a perspective view of three out of the six prism components of the annular expander shown in FIG. 18. As shown in FIGS. 19-21, the annular expander 1748 has six prisms 1748 a-1748 f spaced about a centerline of the annular expander 1748 such that one end of the six prisms 1748 a-1748 f collectively defines an annular-like, multi-faceted surface having a size different from the size of the annular-like, multi-faceted surface collectively defined by the other end of the six prisms 1748 a-1748 f. The light enters the end of the prisms 1748 a-1748 f defining the smaller, annular-like, multi-faceted surface and exits the end of the prisms 1748 a-1748 f defining the larger, annular-like, multi-faceted surface, as shown in FIG. 19. FIG. 20 shows the end view of the annular expander 1748 from the larger, annular-like, multi-faceted surface towards the smaller, annular-like, multi-faceted surface into the page. Although not explicitly shown, the annular reducer 1749 would be configured similar to the annular expander 1748 except that light enters the end of the prisms defining the smaller, annular-like, multi-faceted surface and exits the end of the prisms defining the larger, annular-like, multi-faceted surface.

Returning the FIGS. 17 and 18, the three-phase interferometric microscopy system will now described with reference to light propagating through the system. As shown in FIG. 17, energy is radiated from the sample (labeled “S” in FIG. 17) and reflected by mirror 1744 to beam splitter/annular reverser system 1740. The 50:50 beam splitter 1741 of beam splitter/annular reverser system 1740 (shown in FIG. 17) sends a portion of the light to mirror 1742 and a remaining portion of the light to mirror 1743. To separate the light into two beam portions, the mirror 1742 reflects an outer beam portion C towards mirror 1745 while reflecting an inner beam portion beam B to mirror 1744. Mirrors 1745 and 1744 reflect outer beam portion C and inner beam portion B, respectively, to annular expander 1748 and annular reducer 1749, respectively. Annular expander 1748 expands its received light towards mirror 1746, which reflects the expanded light towards beam splitter 1754 (shown in FIG. 17). Similarly, annular reducer 1749 reduces the radius of its received light towards mirror 1747, which reflects the reduced light towards 50:50 beam splitter 1754.

50:50 beam splitter 1754 reflects half of the light received from mirrors 1746 and 1747 to 66:33 beam splitter 1755 and transmits the other half of the light to mirror 1756. 50:50 beam splitter 1754 also reflects half of the light received from mirror 1743 to mirror 1756 and transmits the other half of the light to 66:33 beam splitter 1755. The 66:33 beam splitter 1755 reflects 33% of the received beam towards 50:50 beam splitter 1758 and transmits the remaining 66% of the received beam towards lens 1760. Mirror 1756 reflects the received beam to 50:50 beam splitter 1758 while also adding an extra phase shift. The 50:50 beam splitter 1758 reflects 50% of the received beam from mirror 1756 towards lens 1780 and transmits 50% of the received beam from mirror 1756 towards lens 1770. The 50:50 beam splitter 1758 also reflects 50% of the received beam from the 66:33 beam splitter 1755 towards lens 1770 and transmits 50% of the received beam from the 66:33 beam splitter 1755 towards lens 1780.

Note that the 50:50 beam splitters 1754 and 1758, mirror 1756 and 66:33 beam splitter 1755 can be arranged such that a difference exists in the light path lengths. More specifically, the optical path distance, x, between the 50:50 beam splitters 1754 and 1758 and the 66:33 beam splitter 1755 can be less than the distance, x+δ_(beamsplitter), between the 50:50 beam splitters 1754 and 1758 and the mirror 1756. Assuming, for example, that the optical path difference between the two distances, δ_(beamsplitter), can be selected to be ⅛ of the wavelength, λ, and the wavelength, λ, is 580 nm, the difference between the two distances, δ_(beamsplitter), can be selected at 72 nm.

The beam sent to lens 1760, the beam sent to lens 1770 and the beam sent to lens 1780 each have substantially equal amplitudes. Lenses 1760, 1770 and 1780 focus the respective beams to detectors 1765, 1775 and 1785, respectively. The intensity of the emitted image on each of the three detectors 1765, 1775 and 1785 changes as a function of the phase/z position because the relative intensities among the three detectors 1765, 1775 and 1785 have about a 120-degree phase shift. FIG. 22 shows a graph of the intensities of the radiated light detected at the three detectors shown in FIG. 17. By monitoring the intensities detected by the three detectors 1765, 1775 and 1785, the phase shift from the initial offset between the inner beam portion B and the outer beam portion C can be estimated as described above in connection with FIGS. 1-11.

This phase shift, δ, can be used to calculate the z position of the emitter: Δz=δλ/(n(1−cos(θ)).

In some embodiments, calibration techniques can be used to correct for non-balanced intensity sensitivities of the detectors and phase mixing angles. Although the embodiments shown in FIGS. 17-22 are described with respect to three detectors having approximately 120-degree phase shifts, other embodiments having a different number of detectors are possible. In yet other embodiments, the imaged spot shape can be used to resolve any cyclic or positional ambiguity. Alternatively, the imaged spot shape can be used as additional information for the positional calculation.

FIGS. 23-25 show an example of an annular expander, according to another embodiment of the invention. More specifically, FIG. 23 shows an end view of the annular expander; FIG. 24 shows a cross-sectional view of the annular expander of FIG. 23 along the line XX-XX; FIG. 25 shows a cross-sectional view of the annular expander of FIG. 23 along the line XXI-XXI. As shown in FIGS. 23-25, the annular expander 2300 has six prisms 2310 a-2310 f, six mirrors 2320 a-2320 f and six mirrors 2330 a-2330 f. The prisms 2310 a-2310 f are spaced about a centerline of annular expander 2300 such that one end of the six prisms collectively defines an annular-like surface having a size different from the size of the annular-like surface collectively defined by the other end of the six prisms. The mirrors 2320 a-2320 f are arranged about the centerline of annular expander 2300 and within an interior central portion defined by the six prisms 2310 a-2310 f. The mirrors 2330 a-2330 f are arranged about the centerline of annular expander 2300 and about an outer portion defined by the gaps between adjacent prisms 2310 a-2310 f.

As shown in FIG. 24, each prism 2310 includes an ingress face 2310′ and an egress face 2310″. The ingress face 2310′ is parallel to and has a smaller inner radius and outer radius than egress face 2310″. Each prism 2310 also includes two side walls 2310′″ and 2310″″ angled relative to the ingress face 2310′ and the egress face 2310″ such that the condition for total internal reflection (TIR) is satisfied. Thus, light enters each prism 2310 at the ingress face 2310′, undergoes a reflection due to TIR at side wall 2310′″ and a reflection due to TIR at side wall 2310″ ″, and then exits from the egress face 2310″. This allows the radius of the incident light to be expanded as it exits the annular expander 2300.

As shown in FIG. 25, the mirrors 2320 a-2320 f reflect the light along and near the centerline of annular expander 2300 outwardly towards mirrors 2330 a-2330 f, respectively. Mirrors 2330 a-2330 f reflect the received light between the gaps defined by adjacent prisms 2310 a-2310 f. For example, mirror 2320 a reflects a portion of light near the centerline of annular expander 2300 towards mirror 2330 a, which reflects this lights between the gap defined between prism 2310 a and prisms 2310 b. Collectively, this allows the central portion of the light to be expanded, providing increased throughput.

Although not explicitly shown, an annular reducer can be configured similar to the annular expander 2310 except that light enters the end of the prisms defining the larger, annular-like surface and exits the end of the prisms defining the smaller, annular-like surface. Although the embodiments shown in FIGS. 19-20 and 23-25 shown with specific refractive and/or reflective components, other embodiments and arrangements are possible including, for example, embodiments having diffractive components or any combination thereof.

Some embodiments include a processor and a related processor-readable medium having instructions or computer code thereon for performing various processor-implemented operations. Such processors can be implemented as hardware modules such as embedded microprocessors, microprocessors as part of a computer system, Application-Specific Integrated Circuits (“ASICs”), and Programmable Logic Devices (“PLDs”). Such processors can also be implemented as one or more software modules in programming languages as Java, C++, C, assembly, a hardware description language, or any other suitable programming language. A processor according to some embodiments includes media and computer code (also can be referred to as code) specially designed and constructed for the specific purpose or purposes. Examples of processor-readable media include, but are not limited to: magnetic storage media such as hard disks, floppy disks, and magnetic tape; optical storage media such as Compact Disc/Digital Video Discs (“CD/DVDs”), Compact Disc-Read Only Memories (“CD-ROMs”), and holographic devices; magneto-optical storage media such as floptical disks, and read-only memory (“ROM”) and random-access memory (“RAM”) devices. Examples of computer code include, but are not limited to, micro-code or micro-instructions, machine instructions, such as produced by a compiler, and files containing higher-level instructions that are executed by a computer using an interpreter. For example, an embodiment of the invention may be implemented using Java, C++, or other object-oriented programming language and development tools. Additional examples of computer code include, but are not limited to, control signals, encrypted code, and compressed code.

While certain embodiments have been shown and described above, it will be understood by those skilled in the art that various changes in form and details may be made. For example, variations in illumination techniques are possible and can allow an improvement to the modulo λ/2 periodicity of the intensity data. In addition, the various embodiments described in the Betzig references can be combined with the various embodiments described herein. For example, in one embodiment, top/bottom TIRF illumination can allow for measurements of a sample having a thickness up to one wavelength, λ. In alternative embodiments, coherent re-phasing of short pulses can form an illumination plane. In yet other embodiments, side illumination to form a Gaussian waist at the sample at various heights in the z-coordinate can further extend the range in the z-coordinate direction without resorting to physical sectioning. Thus, it should be understood that the systems described herein can include various combinations and/or sub-combinations of the components and/or features of the different embodiments described. 

1. An apparatus, comprising: an optical system having a beam splitter device and at least three detectors, the optical system configured to produce a plurality of images of an optical source in a first dimension and a second dimension substantially orthogonal to the first dimension based on an interference of an emission of the optical source in a first direction and an emission of the optical source in a second direction, the second direction being different from the first direction, the beam splitter device configured to interfere the emission in the first direction and having a first source phase and the emission in the second direction and having a second source phase to produce at least three output beams, the at least three optical detectors configured to detect optical energy from the at least three output beams of the beam splitter device, each detector from the at least three detectors being uniquely associated with an image from the plurality of images at a time; and a processor configured to calculate a position of the optical source in a third dimension based on the plurality of images, the third dimension being substantially orthogonal to the first dimension and the second dimension; the interference between the emission of the optical source in the first direction and the emission of the optical source in the second direction being caused at least in part by a difference in path length between the emission in the first direction and the emission in the second direction.
 2. The apparatus of claim 1, wherein the processor is configured to calculate the position in the third dimension based on an amplitude difference between a portion of the at least two images corresponding to the optical source.
 3. The apparatus of claim 1 wherein the plurality of images includes a first image and a second image, the first image being associated with a first portion of the emission from the optical source in the first direction and a first portion of the emission from the optical source in the second direction, the second image being associated with a second portion of the emission from the optical source in the first direction and a second portion of the emission from the optical source in the second direction, the first portion of the emission in the first direction being different from the second portion of the emission in the first direction, the first portion of the emission in the second direction being different from the second portion of the emission in the second direction.
 4. The apparatus of claim 3, wherein: the optical system is configured to align the first portion of the emission in the first direction and the second portion of the emission in the second direction with respect to a first detector from the at least three detectors, and the optical system is configured to align the second portion of the emission in the first direction and the first portion of the emission in the second direction with respect to a second detector from the at least three detectors.
 5. The apparatus of claim 1, wherein the optical system is configured to produce the plurality of images such that each image from the plurality of images is based on a phase of the emission in the first direction at the beam splitter and a phase of the emission in the second direction at the beam splitter independent from the phase of the emission in the first direction.
 6. The apparatus of claim 5, wherein the first source phase of the emission in the first direction at the beam splitter and the second source phase of the emission in the second direction at the beam splitter is based on the difference in path length between the emission in the first direction at the beam splitter and the emission in the second direction at the beam splitter.
 7. The apparatus of claim 1, wherein the optical system includes a first optical component configured to define a first optical path between the optical source and the beam splitter device; a second optical component configured to define a second optical path between the optical source and the beam splitter device; and at least two optical detectors configured to detect optical energy from the at least two output beams of the beam splitter device.
 8. The apparatus of claim 7, wherein the optical system includes at least one spectral filter disposed along an optical path between the optical source and the at least two optical detectors, the at least one spectral filter being a bandwidth filter, low-pass filter or high-pass filter.
 9. The apparatus of claim 7, wherein the optical system includes: a multi-axis positioner coupled to at least one of the beam splitter device, the first optical component, the second optical component or at least two optical detectors.
 10. The apparatus of claim 1, wherein: the at least three detectors is a first set of detectors, the plurality of images is a first plurality of images, the optical source is a first optical source having a wavelength, the optical system further includes a second set of detectors that includes at least three detectors and a plurality of dichroic filters, the second set of detectors configured produce a second plurality of images having a wavelength different than the wavelength of the first optical source, each dichroic filter from the plurality of dichroic filters is uniquely associated with a detector from the first set of detectors and a detector from the second set of detectors, and each dichroic filter from the plurality of dichroic filters is configured to separate optical emissions from the first optical source and optical emission from the second optical source.
 11. The apparatus of claim 1, wherein the optical source is a spatially resolvable point source.
 12. The apparatus of claim 1, wherein the processor is configured to calculate the position of the optical source in the third dimension to sub-diffractive accuracy.
 13. The apparatus of claim 1, wherein the processor is configured to determine the position of the optical source in the first and second dimensions simultaneously with calculating the position of the optical source in the third dimension.
 14. The apparatus of claim 1, wherein the processor is configured to calculate the position in the third dimension based on a relative amplitude between a portion of the at least two images corresponding to the optical source.
 15. The apparatus of claim 1, wherein the optical system is configured to produce the plurality of images using widefield detection of the emission of the optical source in the first direction and the emission of the optical source in the second direction.
 16. An apparatus, comprising: an optical system configured to receive from an optical source at a time an emission having a first source phase and an emission having a second source phase different from the first source phase to produce a plurality of images associated with the time; and a processor configured to calculate a three-dimensional position of the optical source at the time to an accuracy less than a diffraction limit defined by the optical system.
 17. The apparatus of claim 16, wherein the optical system is configured to receive the emission having the first source phase from a direction and the emission having the second source phase from a direction different than the direction of the emission having the first source phase.
 18. The apparatus of claim 16, wherein the processor is configured to calculate the three-dimension position of the optical source based on an amplitude difference between a portion of the plurality of images corresponding to the optical source.
 19. A method, comprising: interfering optical energy emitted from an optical source in a plurality of directions to produce a plurality of output beams; substantially simultaneously detecting a plurality of images based on the plurality of output beams, each image from the plurality of images being uniquely based on an output beam from the plurality of output beams; and determining three-dimensional position information of the optical source based on an intensity of each image from the plurality of images.
 20. The method of claim 19, further comprising: receiving the optical energy such that the optical energy from each direction from the plurality of directions has a source phase different than the source phase of optical energy in at least one different direction from the plurality of directions; and each output beam from the plurality of output beams has a phase different from the phase of the remaining output beams from the plurality of beams.
 21. The method of claim 19, further comprising: the determining including calculating the three-dimensional position information based on an intensity of a portion each image from the plurality of images corresponding to an image of the optical source within the plurality of images.
 22. The method of claim 19, wherein the optical source is switchable.
 23. The method of claim 19, wherein the optical energy emitted from the optical source is a sequence of single photons.
 24. The method of claim 19, wherein the optical source is switchable, the method further comprising: activating the switchable optical source before the interfering; and exciting the switchable optical source after the activating and before the interfering such that optical energy emitted from the switchable optical source is emitted in a first direction and a second direction from a plurality of directions.
 25. The method of claim 19, wherein the optical source is a point source or a quasi point source.
 26. The method of claim 19, wherein determining three-dimensional position information of the optical source includes determining three-dimensional position information to sub-diffractive accuracy.
 27. A method, comprising: receiving a first of input optical beam and a second input optical beam, the first input optical beam from a plurality of input optical beams having a source phase different from a source phase of the second input optical beam from the plurality of input optical beams, the first input optical beam including an image plane and the second input optical beam including an image plane; and producing a first output optical beam, a second output optical beam and a third output optical beam from the first input optical beam and the second input optical beam, the first output optical beam including a first portion of the first input optical beam and a first portion of the second input optical beam and a first phase difference between the first portion of the first input optical beam and the first portion of the second input optical beam, the second output optical beam including a second portion of the first input optical beam and a second portion of the second input optical beam and a second phase difference between the second portion of the first input optical beam and the second portion of the second input optical beam, the third output optical beam including a third portion of the first input optical beam and a third portion of the second input optical beam and a third phase difference between the third portion of the first input optical beam and the third portion of the second input optical beam.
 28. The method of claim 27, wherein the first phase difference between the first portion of the first input optical beam and the first portion of the second input optical beam is different from the second phase difference between the second portion of the first input optical beam and the second portion of the second input optical beam, the first phase difference between the first portion of the first input optical beam and the first portion of the second input optical beam is different from the third phase difference between the third portion of the first input optical beam and the third portion of the second input optical beam.
 29. The method of claim 27, further comprising: determining a location of an optical source based on the first phase difference between the first portion of the first input optical beam and the first portion of the second input optical beam, the second phase difference between the second portion of the first input optical beam and the second portion of the second input optical beam and the third phase difference between the third portion of the first input optical beam and the third portion of the second input optical beam.
 30. The method of claim 27, further comprising tuning each of the phase differences by adjusting a path length of at least one of the portions relative to at least one of the other portions.
 31. A method comprising: directing at least two optical beams emitted from an optical point source along at least two distinct directions, each direction being associated with a distinct optical path; dividing each of the optical beams into two or more divided beams including dividing at least one of the optical beams into three or more divided beams; combining at least two or more of the divided beams to form a plurality of combined beams, each of the combined beams including an interference between the two or more divided beams caused at least in part by a difference in source phase between the optical beams due to a difference in the optical path length of the optical beams; separately detecting an intensity of each of the combined beams; determining a position of the optical source in a first orthogonal dimension and a second orthogonal dimension; and determining a position of the optical source in a third orthogonal dimension based on the separately detected intensities of each combined beams.
 32. The method of claim 31, wherein directing the at least two optical beams includes directing each of the at least two optical beams through an opposing lens.
 33. The method of claim 31, wherein directing the at least two optical beams includes directing each of the at least two optical beams through different portions of a lens.
 34. The method of claim 31, wherein separately detecting the intensity of each of the combined beams includes simultaneously detecting the intensity of each of the combined beams.
 35. An apparatus comprising: a beam splitter configured to receive at least two optical beams emitted from an optical point source along two distinct optical paths, divide each of the optical emission beams into two or more divided beams, and combine at least two or more of the divided beams to form a plurality of combined beams, each of the combined beams including an interference between the two or more divided beams caused at least in part by a difference in source phase between the optical beams due to a difference in the optical path length of the optical beams, at least two detectors configured to separately detect an intensity of each of the combined beams output from the beam splitter; and a processor configured to receive the output from the at least two detectors, and to calculate a position of the optical source in first, second, and third orthogonal dimensions to sub-diffractive accuracy.
 36. A method comprising: activating a statistically sparse subset of switchable optical sources in a sample; exciting the activated switchable optical sources such that optical energy is emitted from the activated switchable optical sources in a plurality of directions; interfering the optical energy emitted from the activated switchable optical sources in the plurality of directions to produce a plurality of output beams; and determining three-dimensional position information of the optical sources based on an intensity of each output beam from the plurality of output beams. 